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We investigate the Kaluza-Klein braneworld cosmology from the point of view of observers on the 
brane. We first generalize the Shiromizu-Maeda-Sasaki (SMS) equations to higher dimensions. As 
an application, we study a (4 + n)-dimensional brane with n dimensions compactified on the brane, 
in a (5 + n)-dimensional bulk. By assuming that the size of the internal space is static, that the 
bulk energy-momentum tensor can be neglected, we determine the effect of the bulk geometry on 
the Kaluza-Klein braneworld. Then we derive the effective Friedmann equation on the brane. It 
turns out that the Friedmann equation explicitly depends on the equation of state, in contrast to 
the braneworld in a 5-dimensional bulk spacetime. In particular, in a radiation-dominated era, the 
effective Newton constant depends on the scale factor logarithmically. If we include a pressureless 
matter on the brane, this dependence disappears after the radiation-matter equality. This may be 
interpreted as stabilization of the Newton constant by the matter on the brane. Our findings imply 
that the Kaluza-Klein braneworld cosmology is quite different from the conventional Kaluza-Klein 
cosmology even at low energy. 

PACS numbers: 98.80.Cq, 98.80.Hw 

I. INTRODUCTION 

Early universe models are usually motivated by theories trying to describe fundamental physics at high energies. 
String theory is a leading candidate to describe high energy physics and since it requires 10 dimensions to be consistent, 
it is natural to consider higher dimensional universes. To reconcile a higher dimensional universe with our empirically 
4-dimcnsional universe, the traditional approach has been to resort to the Kaluza-Klein compactification. Another 
option, the braneworld scenario, was proposed about a quarter of a century ago [j], and has been revived recently, 
stimulated by the discovery of D-branes @. In particular, Randall and Sundrum (RS) proposed an interesting 
framework [j, partly inspired by the Horava-Witten model [||. RS formulated the universe as a domain wall 
in 5-dimensional Anti deSitter (AdS) spacetime. Much work has been done on its cosmology (see @, H, S f° r 
reviews ) and black hole physics [H Gj, 0, HI US Ell • 

However, the RS framework with codimcnsion one braneworld 
is insufficient to reach 10 dimensions. To go beyond five dimensions while keeping our spacetime dimensions to four, 
we need to consider either higher codimensions [TtI . [l8j or Kaluza-Klein compactification on the brane. The former 
option, i.e. to realize a higher codimcnsion braneworld is difficult due to the strong self-gravity of the brane. In fact, 
a higher codimension braneworld develops a severe singularity except for codimcnsion two models. As a result, no 
successful cosmological model is known. Even in the case of codimension two models, it seems almost impossible to 
construct a consistent cosmological model due to the subtlety of the conical singularity 0, [H, [5l|, H3] ■ The latter 
option, i.e. to consider a Kaluza-Klein cosmology on the brane [23|, (2f| [26| , is our concern in this paper. 

One might think it is a trivial task to construct braneworld models with Kaluza-Klein compactification. Unfortu- 
nately, it is not so [27j . In the case of the RS model, the bulk geometry is given and static. Hence, the cosmology on 
the brane is simply due to its motion in the bulk spacetime. In the case of Kaluza-Klein brancworlds, however, the 
bulk geometry is not known a priori (28l. l29l. |30| . Moreover, as we require the internal space to be static, we might have 
to take into account the matter in the bulk, in the form for instance of fluxes. It makes it difficult to solve the bulk 
geometry in most cases. In general, we have to solve the bulk geometry and the brane motion at the same time and 
explicit analytical examples are difficult to construct (see e.g. [3l| for anistropic 5D bulk-brane configurations, with a 
problematics similar in spirit to Kaluza-Klein braneworlds). Although numerical methods seem to be inescapable in 
general, an analytical approach would be useful even if it is a modest one. 

In this paper, we make a first step in the analytical description of Kaluza-Klein brancworlds. Here, we do not 
intend to solve the bulk geometry. Instead, we use the Shiromizu-Maeda-Sasaki (SMS) equation [32| to analyze the 
Kaluza-Klein cosmology. Of course, this effective equation cannot be solved without knowing the projected Weyl 
tensor. Hence, we take the following strategy. We use the staticity of the internal space as a principle to constrain 
the unknown bulk geometry. We also assume that the bulk matter can be neglected, at least in the vicinity of the 
branes, in the regimes which we study. Then, we can determine the Friedman equation on the brane. Interestingly, 
the resultant Friedman equation is found to depend on the equation of state of the matter explicitly. In particular, the 
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effective Newton constant varies logarithmicaily at a radiation-dominated stage. Thus the Kaluza-Klcin braneworld 
cosmology appears to be quite different from the conventional Kaluza-Klein cosmology even at low energy. 

The organization of this paper is as follows. In section 2, we consider a braneworld model with the bulk matter in 
general dimensions and derive the effective SMS equations on the brane. In section 3, we apply the SMS equations to 
a (5 + n)-dimcnsional Kaluza-Klein braneworld model and derive the effective Friedmann equation on the brane by 
imposing the stability of the internal space. The final section is devoted to conclusion. 



II. SMS EFFECTIVE EQUATION IN (D+l)-DIMENSIONS 

In this section, we derive the effective gravitational equations on the brane for any dimension. To be as general as 
possible, we also include a bulk energy momentum tensor. 
The action we consider is 

S=^ J d d+1 x^g[R-2k]-o j d d xJ—g + S m ; A = -^^±± (1) 

where k 2 , i and a are the gravitational coupling constant, the scale of the bulk curvature radius and the tension of 
the brane, respectively. We assume a negative cosmological constant in the bulk. Here, S m represents the action for 
the matter both in the bulk and on the brane. The (d+ l)-dimensional and d-dimensional metrics are represented by 
g and g, respectively. 

We consider a d-dimensional brane with (d — 4) compactificd dimensions. To describe the bulk spacctime, we can 
use Gaussian Normal coordinates so that the metric takes the form 

ds 2 = dy 2 + g^(y, x^dx^dx" , (2) 

and the brane position is y = in this coordinate system. One can deduce the effective equation on the brane following 
SMS. The extrinsic curvature is defined as 

Id 1 

Using the extrinsic curvature, we can write down the Einstein equations in (d + l)-dimensions as 

{d+1) G\ = -\R + \k 2 \K^K a , = + K 2 T% , (4) 

(d+1) G% = -VaV + V M tf = k 2 T\ , (5) 
(d+1 >G% = G\ + (K» u - 5£K) tV - KK% + ^ (K 2 + K a ?K a0 ) = ^"^ ff + K 2 S" v 5(y) + n 2 T\ , (6) 

where G^ v is the <i-dimcnsional Einstein tensor, and T^ v , T ylx , and T yy are the components of the bulk energy 
momentum tensor. Here V M denotes the covariant derivative with respect to the metric g^, and S^ u = —crg^ + 
is the energy momentum tensor on the brane, where is the energy momentum tensor of the brane matter other 
than the tension. Then, the junction conditions are given by 

[K» v - 5£K] \ y=Q = y i-aSH + P„) , (7) 
where we have assumed ^-symmetry. Combining Eqs. (j4| with ([6|), we have 



d- 1 



d(d-l) 

The trace of this equation gives 

d nd-2 



^g,uK 2 + l -g^K^K aP - - J^-j (t^ - ^i^lO • (8) 



K = K af3 K r- — + k 2 - -T y —T» (9) 
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Also the following components of the Weyl tensor are relevant. 

1 f u 1 \ d-2, d-2 d-'^ r- x K 

Lw ~ d-i \ ^ d 9flv J + d - 1 ^ y d(d - i) 9 ^ K -y + d - 1 M Al/ 

+ - d 9^K a0 - • (10) 

The above components of the Weyl tensor can be rewritten by using Eqs. ([8]) and (|9j) as 

W d -^g, v T\ + ^Vt»„) . (11) 

Thus off the brane, using these components of the Weyl tensor, Eq. © is expressed as 

= —C yfiyt , — K^Kxv + KK^ + -g^K al3 K a p — -g^K 2 + ^j- -<? M „ 

+ jk 2 - ^g^T a a ^j + ^-^ K 2 T y v9tlv (12) 

where we stress that the term T a a is the trace defined with respect to the d-dimensional metric g, and not the full 
trace defined with respect to g. Eliminating the extrinsic curvature by using the junction conditions ijTj). and assuming 
the RS type relation 

, 2(d-l) 

k a = , (13) 

we finally obtain the <i-dimcnsional generalization of the SMS equations, 

R»v - \g^R = -E^ + 8ttG tfj, v + + ^— jK 2 (t^ - \g^T a ^}j + d —^L K ^ v g tiV , (14) 

where we have defined the Newton constant in (^-dimensions by 

(d-2)n 2 , . 

= 1 J , (15) 
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the tensor 

1 1. >. 1 



- 4(^ _ i} U ^ u ~ 4^ A *A!/ + -9^t a ^t a i3 - - ^g^t 2 , (16) 



and the projected Weyl tensor 



E(j, U — Gy^yv | Q . (17) 



These results are also obtained in [33j, l34[ . 

From the Bianchi identity satisfied by the Einstein tensor, we can deduce a constraint equation on the tensors that 
appear on the right hand side of 

- k 4 V%„ + ^jVVT^ - J0^ K2 V„T a a + 1_^ K 2 V„T% , (18) 

where we have assumed the conservation of the energy-momentum tensor for the matter on the brane 

V% v = , (19) 

i.e. we forbid the possibility of energy exchange between the brane and the bulk (as studied e.g. in [35| in 5D brane 
cosmology). We also have the conservation law for the bulk matter, which can be decomposed as 

= d y T y y - KT\ + K\T\ + V M T"„ , (20) 
= a y T% - KT\ + V„T% . (21) 
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Of course, the above equations do not form a closed system, because we do not know E^ v . In other words, without 
knowing the bulk geometry, we cannot solve the effective Einstein equations (|14[) . However, we may regard the SMS 
equation as an initial value equation. Namely, once is given from the SMS equation, we can solve the (d + 1)- 
dimensional Einstein equations along the y-direction to obtain the bulk geometry. In this picture, unless we impose 
some conditions on the properties of the spacetime, there will be too many allowed bulk solutions, and most of the 
solutions will be physically meaningless. 

In the next section, we will try to solve the effective Einstein's equations to study the cosmology of Kaluza-Klcin 
brancworlds from the SMS equation, by assuming that the internal dimensions are static and that the bulk energy- 
momentum tensor can be neglected on the brane. Given these conditions, we will show that, similarly to the cosmology 
of a codimension 1 brane in an empty 5D bulk [36|, [37j, the effective Friedmann equations can be solved, up to a 
constant of integration. However, in contrast to the 5D bulk, we will not give a bulk geometry associated with this 
cosmology. 



III. KALUZA-KLEIN BRANEWORLD COSMOLOGY 



To realize a braneworld in higher dimensions, it seems natural to consider a codimension one brane with internal 
dimensions that are compactified a la Kaluza-Klein, in short a Kaluza-Klein braneworld. Here, for simplicity, we 
consider n internal dimensions compactified on a torus. The brane thus represents a (4 + n)-dimcnsional spacetime 
embedded in a (5 + n)-dimcnsional spacetime. Since we wish to study the cosmology, we consider metrics of the form 



ds 2 = -dt 2 + a 2 (t)S tj dx i dx j + b 2 {t)5 a pdz a dz l: 



(22) 



where the x l are the three ordinary spatial coordinates and the z a are the internal coordinates. For simplicity, we 
assume that there is a single scale factor b characterizing the size of the internal dimensions. The scale factor a is the 
usual scale factor for the external space. 

We can imagine two kinds of matter, the matter in the bulk and on the brane. The bulk matter is important to 
get a well-behaved geometry in the bulk [39L |40| . However, for simplicity, we suppose here that we can ignore it for 
the cosmology on the brane. Hereafter, we will consider only the matter on the brane. Note a recent work, where 
a similar analysis was considered for a 6D Kaluza-Klein brane embedded in a 7D bulk spacetime, and which takes 
into account a bulk energy- momentum tensor and the possibility of brane-bulk energy exchange [381 ] . Because of the 
symmetries, the energy-momentum tensor is restricted to be of the following form, 



tpv = (p, Pa 2 5ij,Qb 2 5 a 



(23) 



where p is the energy density, P the external pressure and Q the internal pressure. Similarly, the projected Weyl 
tensor is of the form 

E„ v = (e, X a 2 S. u ^b 2 S a p) . (24) 
+ 3% + n£ = 0. The components of the quadratic 



Moreover, the traceless property of E^ v implies the relation 
tensor in the energy-momentum tensor p^|) are given by 



1 



TOO 



TTq/3 



8(3 + n) 
1 

8(3 + n) 
1 



(n + 2)p 2 - 3n (P - Q) s 

[(n + 2)p 2 + 2p (2P + nQ) + n (P - Q) (P - 3Q)] a 2 8 t 
[(n + 2)p 2 + 6pP + 2{n - l)pQ + 3 (n(P - Q) 2 - 2Q 2 + 2PQ)] b 2 5 a(3 . 



(25) 
(26) 
(27) 



8(3 + n) 

Substituting the metric (|22|) and the tensors (|23|). (f24|) and (|25H27j) in the effective Einstein equations (fT4|) . one finds 

n(n-l) „ 2 



mt + 3nH a H b 



-2H a - 3H 2 - 2nH a H b - nH b - 



-Hi = 8nGp 



n(n + 1) u2 



8(3 + n) 

^4 



Hi = 8nGP 



(n + 2)p 2 - 3n (P - Q f - e (28) 



[(n + 2)p 2 



-3H a - m 2 a - 3(n - l)H a H b - (n - l)H b 



ti h 



8(3 + n) 

2p(2P + nQ) + n(P -Q)(P -3Q)]-x , (29) 

„4 



8itGQ 



8(3 + n) 



[(n + 2)p 2 + 6pP + 2(n - l)pQ 
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+3 (n{P - Q) 2 - 2Q 2 + 2PQ)] -£ . (30) 
In addition to the above equations, we need the conservation law for the matter (|19p . which is given here by 

p + 3H a (p + P)+nH b (p + Q) = . (31) 
The constraint equation for the projected Weyl tensor (TT8|) in the absence of the bulk matter can be written as 

3n r • 1 

k + 3H a (e + x)+nH b {e + t) + ——K 4 (P-Q) P - Q + H a (p + P) - H b (p + Q) = 0, (32) 

4(3 + n) L J 

where we have defined the Hubble parameters H a = a/a and H b = b/b. In order to integrate explicitly these equations, 
we will assume simple equations of state for the anisotropic fluid, namely P = wp and Q = vp, where w and v are 
constants. Then the equations (|28ti30|) become 

3H 2 a + 3nH a H b + ^""^ ffg = SnGp + * + - {2 + n(l - 3{v - w) 2 )} p 2 - e ,(33) 

-2H a - 3H 2 a - 2nH a H b - nH b - " ( " + 1} H 2 = 8nGwp + ^ - {2 + Aw+ 

2 8(3 + n) 

+n(l + 2v + 3v 2 - Avw + w 2 )} p 2 - x , (34) 

-3H a - 6H 2 - 3(n - l)H a H b - (n - l)H b - ~ ^ H? = 8t:Gv P + {2 - 2v - 6v 2 + 6w + 6vw+ 

2 8(3 + n) 

+n(l + 2v + 3v 2 - 6vw + 3w 2 )} p 2 - - (e - 3*) (35) 

' n 

and the conservation law reduces to 

p + 3(l + w)H a p + n(l + v)H b p= . (36) 

The constraint (13"2l is written as 



e + 3H a (e + X ) + H b {{n + l)e - 3 X ) + .J" [(w - vfpp + (H a (l + w) - H b (l + v))(w - v)p 2 ] = .(37) 

4(3 + n) 

Since we do not know e, \ an d £> we cannot solve the above equations. We need to solve the bulk geometry to 
obtain e, x an d £, in general. Here, instead of solving the bulk geometry, we impose the stability of the internal space, 
that is, we put 6=1. By doing so, we may have a singularity in the bulk. However, if we allow the existence of 
matter in the bulk, it is reasonable to expect that the bulk geometry can be made regular by a suitable choice of the 
bulk matter. What kind of matter is necessary is a different issue. Here we assume the staticity of the internal space, 
and seek for an effective Friedman equation on the brane. Under this staticity assumption, Eqs. (|3"3"|) - (|3"T|) reduce to 

3H 2 = SirGp + f {2 + n{\ - 3(v - w) 2 )} p 2 - e , (38) 
8(3 + n) 

-2H a - 3H 2 = 8irGwp + — {2 + Aw + n(l + 2v + 3v 2 - Avw + w 2 )} p 2 - x , (39) 

8(3 + n) 1 ' 

K 4 1 

-3H a - m 2 = 8irGvp + — {2 - 2v - 6v 2 + 6w + 6vw + n(l + 2v + 3(v - w) 2 )} p 2 - - (e - 3 X ) , (40) 

8(3 + n) 1 n 

and 

p + 3(l + w)H a p = , (41) 
3n 

e + 3H a (e + X ) + — - k 4 [(w - vfpp + H a (l + w)(w - v)p 2 ] = . (42) 
4(3 + n) 

What we want is the effective Friedman equation in the Kaluza-Klcin brancworld. For that purpose, a glance of 
Eq. (|38[) tells us that we need to know e, which is a component of the projected Wely tensor which encodes some 
information about the bulk geometry. For general w, Eq. (|4Tj) is solved to give 

p = a- 3 ^™) , (43) 
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where we have absorbed a constant factor into the scale factor by rescaling it. This is a standard result. Eliminating 
H and H 2 from Eqs. (|3"5 j) -([l0" ]) . we obtain 

- n£ = 3% - e = 8ttG-^— (3w -2v- l)p + n 4 —^ — -v (1 - 3w + 3v) p 2 . (44) 
s 2 + n y 4(3 + n) v ' y ' 

Thus the components of the projected Weyl tensor are related to the matter on the branc. Substitution of the above 
result (|44p into the constraint equation for the dark radiation (|42|) gives the equation 

e + AH a e = -8nG-?—(3w - 2v - l)H a p - k 4 A . - , (1 - 3w + 3v) {w + 3(1 + - «)} H a p 2 . (45) 
2 + ?i 4(3 + n) 

This equation can be integrated easily. Note that the cases w = 1/3 and w = —1/3 need to be considered separately. 
First, we consider the generic case w ^ 1/3, —1/3. The solution is given by 

e = _gg + 87rg n 1 + 2^-3^ q _ 3(1+M)) 
a 4 2 + n 1 - 3u> 



4 



(1 - 3w + 3v) 



8(3 + n) (1 + 3w) 



{w + 3(1 + w){w - v)} a~ 6{1+w ^ , (46) 



where C is a constant of integration which can be interpreted as dark radiation [37], EH . Substituting Eq. (j46|) into 
Eq. ([55|) . we obtain the effective Friedman equation 



where the coefficients are given by 



2{1 -3w-nv) 

Goff = (2 + n)(l-3«,) G ' (48) 

A = 2 + 6^-f n(l + 2r, + 3(r,-^) 4 

24(3 + n)(l + 3w) V ; 

The above equations include the well-known 5D case, corresponding to n = and for which G c s = G and A = k 4 /36. 
By contrast, in higher dimensions, the effective Newton constant that we have defined depends on the equation of 
state. It means that the Kaluza-Klein braneworld cosmology, provided our assumptions are valid, is different from the 
conventional Kaluza-Klein cosmology even at low energies. For w = 0, we have to assume nv < 1 in order to have the 
positive effective Newton constant. It should be noticed that the effective Newton constant becomes negative in the 
regime < nv < 1 and (1 — nv)/3 < w < 1/3. This indicates some transient instability around the matter-radiation 
equality. 

Before proceeding, however, it is important to note that in addition to the assumption of staticity of the internal 
space, we have also assumed that the bulk matter can be neglected and that there is no explicit coupling between the 
matter on the brane and the matter in the bulk. If we relax one or several of these assumptions, the conclusion will 
be altered. 

Now, we consider the radiation-dominated case w = 1/3. In this case, the solution reads 



3C 2tl logo 4 

— r + 87rG v — -. K 

a 4 2 + Ti a 4 16(3 + n) 



irG- v-^f - k 4 — -v(9v ~ A)a~ 8 . (50) 



There appears a logarithmic correction in the above. Substituting Eq. (f5D|) into Eq. (|3"5|) . we obtain the effective 
Friedman equation (|4T[) with the coefficients given by 

Geft = G(l-^vlog-) , (51) 
\ 2 + n a* / 

_ 12 + n(4 + 9^) 

A 144(3 + n) ' (5 } 

where a* is a constant of integration corresponding to the dark radiation component C. Remarkably, the effective 
Newton constant depends logarithmically on the scale factor. This interesting behavior of the cosmological evolution 
occurs only during a radiation-dominated stage. Furthermore, depending on the value a*, the effective Newton 
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constant may become negative. This implies that the dark radiation component should be chosen appropriately in 
order to realize a sensible cosmology on the brane. It may be mentioned that this behavior is in contrast to the case of 
the Brans-Dickc cosmology in which the effective Newton constant begins to depend on time after pressureless matter 
starts to dominate. In a sense, one may say that dark matter stabilizes the effective gravitational coupling constant. 
The relevant constraint comes from nucleosynthesis [42|, [43| : 

AG c ff „ +0.20 



G, 



off 



O.OlloX (53) 



at one-sigma confidence level. It is easy to see that this constraint is satisfied for a sufficiently large a*. 
Finally, we consider the case w = —1/3. The solution is given by 

3G „ ri , , a n , . , log a 

e = -^ + 87rf W X + ^ + K 12(37^ + ^ ■ (54) 



where a logarithmic term appears again. Substituting Eq. (|54|) into Eq. (|38|) , we obtain the coefficients in the effective 
Friedman equation as 

Goff = %- G , (55) 

2 + n 

a 6 + n(2 — 6v — 9v 2 ) A n , a . 

A = - s 2 + 3v 2 log - , 56 

72(3 + n) 36(3 + n) a* 

where a* is the constant of integration we mentioned before. In this case, we have a logarithmic scale factor dependence 
in the coefficient of p 2 . This can have some impact at high energy. Therefore, its effect on the inflationary scenario 
may be interesting. 

Apparently, the Kaluza-Klcin cosmology we have obtained is different from the conventional Kaluza-Klein cosmol- 
ogy. In particular, the effective Friedman equation depends on the equation of state of the matter explicitly. This 
result will hold in more general Kaluza-Klein spacetimes. The reason is the following. The constraint equation (Tig]) 
in the absence of the bulk matter reads 

WE^ = f^W^ . (57) 

Hence, in general, the projected Weyl tensor is affected by the energy momentum tensor on the brane. If the brane was 
isotropic and homogeneous, the matter part would have the additional property that V7r M „ = 0. In our example, this 
can be seen by setting P = Q in Eq. (|3"2"j) . The effect of matter would then not appear in Eq. Ij5"7"|) . This is related to 
Birkoff's theorem. Because of spherical-like symmetry, one does not see the details of the matter contents, but rather 
see the dark radiation which is independently conserved. On the contrary, in the case of Kaluza-Klein braneworlds 
this no longer happens because of the anisotropy of the brane. In particular, ir^ is no longer conserved, =/= 0. 

Thus the anisotropic brane will deform the bulk geometry nontrivially. Then the back-reaction of this to the brane 
will modify the effective Friedman equation. Because of the strong anisotropy of the Kaluza-Klcin braneworld, we 
then expect that the modification of the Friedmann equation will persist even at low energies. 

Since the SMS equation is not a closed system, we cannot formulate cosmological perturbation theory without 
further information. As we noted before, we may regard the SMS equations as giving the "initial data" to solve the 
bulk geometry. Of course, in general, there would be a singularity in the bulk. However, if one introduces some bulk 
matter, one may obtain a non-singular geometry in the bulk. With some explicit bulk matter, one can also formulate 
explicitly the cosmological perturbation theory for the corresponding Kaluza-Klein braneworld. Since there is the 
possibility of a deviation from Newton's law at low energy, it would be interesting to perform this program explicitly 



IV. CONCLUSION 



We have considered a (d + l)-dimensional gravitational system with bulk matter and a <i-dimensional brane, and 
derived the effective d-dimensional Einstein equations on the brane. As an application, we have studied the cosmology 
of Kaluza-Klcin braneworld, with n internal toroidal dimensions. By neglecting the bulk matter and imposing the 
staticity of the internal space, we have obtained a closed set of equations, from which we have been able to derive 
the effective Friedmann equation. We have found that the resultant Friedmann equation explicitly depends on the 
equation of state. In particular, if the universe is dominated by radiation, a resonant contribution to the projected 
Weyl tensor gives a time variation to the effective gravitational coupling constant. This time dependence disappears 
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after the radiation-matter equality, which can be interpreted as a stabilization of the Newton constant by the matter 
on the brane. It should be emphasized that the Kaluza-Klein braneworld cosmology is quite different from the 
conventional Kaluza-Klein cosmology even at low energy. It is in contrast to the fact that the conventional Fricdmann 
equation can be recovered at low e nerg y in the RS braneworld cosmology. Hence, it is important to see whether 
Newtonian gravity can be recovered [44J on the Kaluza-Klein braneworld. This question has already been studied in 
some specific 6D models based on flux compactifications [13, Ef| . 

We have also discussed the method to obtain the bulk geometry from the brane data. It seems always possible to 
adjust E^u so that the braneworld has a static internal space during the cosmic expansion. Of course, in general, the 
geometry of the bulk will be contrived and it will perhaps contain singularities. However, assuming the presence of 
matter in the bulk, there is a chance to have a non-singular bulk. Admittedly, it is a non-trivial problem to find an 
appropriate matter which can stabilize the braneworld without any pathology in the bulk. An alternative approach 
is to start from known bulk solutions in which a brane is embedded. This method has been used very recently to 
study the cosmology of Kaluza-Klein branes, with one internal dimension, in 6D bulk solutions of Einstein-Maxwell 
equations [H,[!3]. 

There are many other issues to be explored. It is interesting to formulate the quantum creation of the Kaluza-Klein 
braneworld as is done in RS models 1481. |49| . It is also important to understand the low e nerg y description of the 
Kaluza-Klein braneworld [lO, HH, l52l. l53l. lEH, [H, l56l and the Kaluza-Klein corrections [13, l58| . It is intriguing to 
consider born-again scenario in higher dimensions [59l. [60| . We leave these issues for future studies. 
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